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Abstract 



Interaction of a fast shock wave generated during a supernova explo- 
sion with a magnetized star-companion of the supernova precursor pro- 
O ' duces a current sheet. We consider the evolution of this current sheet 

and show that a singularity (shock) is formed in finite time within the 



C/3 . ideal MHD framework. Charged particles (electrons) are accelerated in 

^ I ' the vicinity of the singularity, and their distribution function has a plateau 

up to the energies of the order of IQ'^mc? . These fast particles radiate in 
the 7-range in the strong magnetic field of the current sheet {B ~ 10*' G). 
^ ' Radiation is concentrated within a narrow angle around the current sheet, 

, A6' ~ 3 • 10~*, and its spectrum has the maximum at several hundreds 

\l ' of keV . Presented calculations confirm the model of cosmological GRBs 

(N ■ 

proposed by Istomin & Komberg (2002). 

O ■ 1 Introduction 

oo ■ 

' Observations show that some gamma-ray bursts (GRBs) occur at cosmologi- 

cal distances, and some of the bursts are related to the supernova explosions. 
)«j ■ If the radiation of a cosmological GRB is isotropic, its energy release is esti- 

. mated to be 10^^"^'' erg at the hundreds of keV range. Massive supernovae 

' emit comparable amounts of energy, mainly in the form of neutrinos. Thus a 

question about the origin of GRBs' energy supply arises. Connection between 
the GRB phenomenon and that of a supernova seems to be supported by several 
facts. First, the paper by [S] reports the temporal and spatial concurrence of 
the "compact" supernova SN 1998 bw (10*^ erg, type SNIc) with GRB 980425, 
identified with a nearby galaxy which has Z — 0.0085. Secondly, "hills" are 
observed in the light curves of some optical transients several days after the 
GRB; they can be interpreted as a contribution from the supernova burst si- 
multaneous with the GRB. Thirdly, narrow X-ray lines of ions characteristic 
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of supernova shells are detected in the spectrum of the afterglow [S], with the 
velocity of the line-emitting material reaching 0.1c. Obviously, the "energy cri- 
sis" arises if we attempt to relate an isotropic GRB to a supernova. One way 
to avoid this crisis is to suggest that GRBs are associated with strongly coUi- 
mated relativistic ejections radiating within 7-range. This makes it possible to 
decrease by many orders the luminosity from cosmological GRBs corresponding 
to the observed fluxes. [T] proposed this connection between GRBs and strongly 
anisotropic supernova bursts. Summarizing, we see that current understanding 
of the nature of GRBs requires a combination of two conditions: a supernova 
explosion (possibly of special type) and the formation of a narrow (opening 
angle less than a degree) beam of relativistic particles radiating within X- to 
7-range. [5] proposed a model where a beam of relativistic particles is formed 
by the interaction of a fast shock wave generated during a supernova explosion 
with a magnetized star-companion (a neutron star or a white dwarf) of the 
supernova precursor. The flow of the shock matter around the magnetosphere 
of the star-companion forms an almost parallel magnetospheric tail. Magnetic 
field in the tail becomes super strong (about 10® G) due to the cumulative effect. 
Virtually any star-companion interacting with a supernova shock will produce 
such strong magnetic field, due to amplification during the explosive compres- 
sion of conducting medium - the mechanism proposed by to produce very 
strong magnetic fields. For the typical shock parameters, the density psh is of 
order lO^^g/cm'^, and the velocity Vsh is about 4 • lO^cm/ sec, the value of the 
amplified magnetic field is B = {^-k pshV^j^Y^"^ ~ 10® Gauss. For a usual star 
like the sun (which has magnetic field of the order of lOG), the compression 
will be about 3 • 10^ times and the transversal size of the tail will be equal to 
Rq/?> • 10^ ~ 10* cm. Hence, virtually any explosion of a supernova which had a 
star-companion, will result in a highly magnetized elongated plasma formation 
in the shadow of the star-companion. The magnetic field changes its direction 
from one edge of the tail to the other, which leads to the magnetic reconnection, 
magnetic energy release, and acceleration of charged particles. Similar processes 
take place in the magnetospheric tail of the Earth formed by the flow of the so- 
lar wind. In a strong magnetic field these relativistically accelerated particles 
produce synchrotron radiation concentrated in a narrow cone ~ 3 • 10~^. 
The frequency of supernova explosions has the order of 3 • 10^^ supernovaes/yr 
per a galaxy, which gives (3 • 10"'*)^ • 3 • 10~^ • 10^^ ~ 300 radiation bursts per 
year, directed towards the Earth (the Universe has ~ 10^^ galaxies). All this 
may provide an insight into the nature of GRBs. In this paper we provide an 
accurate quantitative theory of the outlined phenomena, namely the evolution 
of the current sheet, the acceleration of charged particles, and the formation of 
the radiation spectrum. 

2 Stationary current sheet 

Magnetic field of a star is almost dipolar at the large distances from the center 
of the star. Its direction changes from one side of the sheet to the other (with 
respect to the axis of the tail) . Interaction of a dense plasma fiow with this dipole 
magnetic field produces a long current sheet. A similar magnetospheric tail is 
formed by solar wind in the dipolar magnetic field of the Earth. Such situation 
may also arise in a close binary system of stars, consisting of a magnetized star 
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Figure 1: Stationary current sheet 

and a massive supernova precursor. A shock wave from the exploding massive 
star interacts with the magnetized star and forms a current sheet, resuhing in 
a GRB. Such model of GRBs was proposed in the papers of [5]. Yet another 
example of the current sheet formation comes from a binary system of neutron 
stars, one of which is a radio pulsar. Relativistic wind flow of the radio pulsar 
around the other star forms a magnetic tail. Discovery of a close binary system 
of neutron stars containing two radio pulsars PSR J0737-3039 A,B was reported 
in the paper of [iTj . Observations of this system show that the wind flow of 
the more powerful pulsar 'A' around the magnetosphere of the radio pulsar 'B' 
takes place inside the light cylinder of the pulsar 'B'. Thus, a long and almost 
parallel tail of magnetic field is formed in the magnetosphere of 'B'. 

To study the phenomena taking place in the magnetosphere of a magne- 
tized star, we assume the sheet to be flat and infinitely long. Therefore, all 
characteristics depend only on one coordinate z, which is perpendicular to the 
current sheet. Our approximation is one-dimensional, which makes possible the 
solution of non-stationary problems. The geometric structure of this problem is 
presented in Figure 1. 

Non-uniform magnetic field is parallel to the x axis and the tail, the y axis 
is perpendicular to the x-z plane. Magnetic field depends linearly on the z 
coordinate 

B = (A'z,0,0), 

where iiT is a constant determined by electric currents flowing in the y direction. 
A particle with charge q and mass m in such non-uniform magnetic field is 
subject to the gradient drift ^4] 



2mqKz^ 



Here v± and p± are the velocity and momentum components of the particle 
perpendicular to the magnetic field. Equation ([1]) is not applicable near the 
plane z ~ 0. The Larmor radius of the particle is greater than the inverse 
gradient of the field there, so the drift approximation is not valid when Kz^ < 
cp±/q. Electric current jd produced by the drift motion of charged particles 
with density n is given by 
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and is independent of their charge. Brackets < ... > mean averaging over the 
distribution function of the particles. In addition to the drift current resulting 
from the drift of the particles in the non-uniform magnetic field, there also exists 
a magnetization current jm due to non-uniform distribution of the magnetic 
moments of the particles, jm — c ■ curZM, M = n</i>, /i = ~'B{p±vj_/2B'^). 
The sum of the drift current and the magnetization current depends only on 
the gradient of the transversal energy density of particles 

c d f < pj_vj_ >\ 

= -bTz r — 2 — ) ■ 

The total current j is formed by the electron and ion flows. Assuming plasma 
to be quasineutral and consisting of electrons and positively charged ions with 
mean charge Ze, we see that the electric current is 



c d 

^y^~Bd^ 
On the other hand, we have 



n.il + Z)^J^-. (4) 



c cK 
—curia = - — 



Finally, 



f- + n,(l + Z) < const = Po. (5) 

Equation (0) is virtually identical to the equilibrium condition of an inhomoge- 
neous plasma in a non-uniform magnetic field, when the equilibrium is supported 
by external pressure Pq. 



3 Evolution of the current sheet 

Now we add an extra assumption that the current sheet is compressed by some 
external force acting along the z axis, and consider the evolution of such sheet. 
We may assume that the evolution is fast enough and it can be described by 
ideal magnetohydrodynamics (MHD). MHD equations for the one-dimensional 
case (all parameters depend only on the z coordinate) are 

dB dB dv^ ^ 
at oz oz 

dt ^ ^ dz ~ ^' 

(dv. dvA^dP I dB 

dp dp dv, 

+Vz-^ + P-^ = 0, 
dt dz dz 

dP dP ^dv, „ 



dt dz d 



z 



During the one-dimensional compression (or expansion) of the current sheet 
magnetic field remains one-dimensional, 

B = (B(z,i),0,0). 
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In addition to the y component of the velocity v, the z component becomes 
time-dependent as well 



V = {Q,Vy{z,t),V:,{z,t)) . 

The plasma density is p = p{z,t), its pressure is P = P{z,t). We use the 
equation of state P = C2P^ , where F is the adiabatic exponent. From the first 
and the fourth equations of system ([6]), we obtain 

dB dB dvz dp dp 

which leads to the following condition 

dB dp 
P^=^di' 

which means that the magnetic field is "frozen" into the plasma. Hence, we have 
B (X p, so let B — Cip. Now, the third and the fourth equations of system ^ 
describe the evolution of the velocity Vz{z, t) and plasma density p(z, t) 

dvz , 9vz\ ^ „ ^ T-idp , Clpdp 
at az J oz 47r oz 

dp dp dvz , , 

^+w^^+p^ = 0. 7) 
ot oz oz 

The change of variables {p, Vz) (t, z) reduces the system of non-linear equa- 
tions ^ to the following linear system 

^^y^ + c-r ^-^— + ^—-0 (8) 

dp dp dv Att dv ' 

dz dt dt ^ 
dv dv P dp 

To simplify the notation we drop the index of Vz- Eliminating the derivatives 
of z we obtain one equation for t{p, v) 

PW^^-""'^' dv^-^dv^^'- ^''^ 

Solution of this equation cannot be expressed in terms of known special functions 
for the arbitrary values of F ^ 2. However, we can obtain a quasiclassical 
approximation to the general solution after applying Laplace transform to the 
variable v. When w > 0, the quasiclassical solution to ([TO)) is 

„ 1 1 r fM , ,,,, 

Heree(p) = (C2Fpi^-3 + C2/(4^p))'/'= (FP + SV(4^))'/%-3/2 = c,/p, and 
Cs is the fast magnetic sound speed. The function f(p) appearing in the solution 
depends on the initial conditions. We may assume that the velocity v is zero 
when t — 0. Accordingly, dv/dz\t=o = as well. It is convenient to write the 
initial condition for the density p at < = in the following way f{p)dp = z, 
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where po is the density at z = 0. From the continuity equation of system ([7]) 
it follows that dp/dt\t=o = 0. From the first equation of the same system (the 
equation of motion) we see that dv/dt\t=Q = — p£,^ (p) / f (p) ■ After we change 
variables to {t(v,p), z{v,p)), our initial conditions become 



dt 




f{p) 


dv 




peip)' 


dz 




0- ^ 


dv 


v=0 


dp y 



dp 



= 0; 



ti=0 



= /(p); tl=o = 0- 

11=0 

In the obtained solution, p+ satisfies the equation v — ^{pi)dpi, and p_ 

satisfies the equation v — S,{pi)dpi. For v > Jq S,{pi)dpi the value p_ 

cannot be determined, because at the limiting value v = Jq £,{pi)dpi the density 
P- vanishes and cannot be negative later. Since ^ = Cs/p, the condition v > 
Jq ^{pi)dpi is virtually identical to t; > Cg, and quasiclassical solution (jlip is 
not valid. 

Substituting pip into ^ gives the following approximation for z(v,p) 



z{v,p) = vt{v,p) 



d_ 

dp 



p / t{vi,p)dvi 



f{pi)dpi 



(12) 



The integrals appearing in expressions (|llll2p have their simplest form when 
r = 2, which is close to the adiabatic exponent for an ideal gas F = 5/3. In this 
case the solution is 



t{v,p) 



4/o k 



1 



8/o 



(13) 



Here we have chosen the function f{p) to be the power function f^p^ where k is 
an arbitrary value and a = 2(26*2 + Cj/47r)i/2. Solutions t{v, p) and z{v, p) ^ 
are always single- valued functions of (v,p), but the functions v{t,z) and p{t,z) 
may not be. The quantity p{t, z) can be found by reversing the second equation 
of system (jl3p . This can be uniquely done until the moment when the shock 
occurs and the Jacobian of the transformation of the unknown variables vanishes 



J 



d{t,z) 
d{p,v) 



9 P 



1 



2fc(l-fc) — 
a^p 



0. 



For fc < and fc > 1 the Jacobian J always vanishes at a value q — l{c?p). 
When the shock occurs, the derivatives dp/dz and dv/dz become infinite which 
leads to the singularity in the electromagnetic field. The question whether 
there is a shock for almost arbitrary initial perturbations is very important, so 
we would like to obtain an exact solution of the MHD equations, and it is only 
possible for F = 2. This solution can be expressed in terms of hypergeometric 
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Figure 2: —?>a'^ J / Af^p^*^ of q at fc = —3/4, the dashed hne is the quasiclassical 
approximation 



functions 



HP>") = 2 



■fp -l-fc,-fc,-,— 

\ 2 a^p / 



'2'a2py /j + l ' 

The hypergeometric function 2Fi{a,b,c,q) appearing in these equations has a 
branch point at q = 1 , which means that solution (I14p cannot be continued into 
the region q > 1. 

Note that quasiclassical solution (fT5|) for F = 2 corresponds to the first two 
terms in the expansion of exact solution ()14p into powers of q. 

Both in the exact solution and in the quasiclassical one, the singularity 
occurs for < g < 1 when the Jacobian J vanishes (the shock point). However, 
in the exact solution the shock occurs only for negative values of k, namely for 
k < —1/2. These values of k correspond to the more natural initial conditions 
when the initial value of the acceleration of the sheet dv/dt\t=o (x does not 
have a singularity at p — Q. The behavior of the Jacobian J oi q = v^/ {a^p) for 
the exact and quasiclassical solutions is shown in Fig. 2. 

The figure shows that the approximate solution describes the behavior of the 
sheet well right up to the shock point, so later we will use only quasiclassical 
solution (fTSll . 

An expression for the y component of the plasma velocity should be added to 
formulas p^ . It can be easily found from the equation dvy/dt + vdvy/dz = 0, 
whose solution is an arbitrary function of the variable p^^^rj, Vy = ^{p^^^rf), 
where 

r,= l+ + 2)«^ p-2fc-42 ^ ^ ^ ^ (15) 

The function $ is determined by an initial condition for the motion of the sheet 
in the y direction at t = 0. The initial velocity Vy can be arbitrary for an 
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infinitely wide sheet. But since the real geometry of our problem is not flat, we 

may assume that Vy ^ = v. 

The electric and magnetic fields are given by 

ft - fl_r n - « _ fe _vy _<i> 

B c 4/oc B c c 

and the other components vanish. 



4 Acceleration of charged particles 

Evolution of the current sheet causes a drastic change of the magnetic field, 
its energy, and the appearance of the electric field. Due to the high speed of 
this process (it lasts for only about L/cs ~ O.lsec), the energy of the magnetic 
field is passed to the energy of the accelerated particles. Acceleration is most 
effective in a region with the fastest evolution, i.e. near the singularity. Let 
us consider acceleration of the charged particles in the vicinity of singularity. 
We assume that these particles are relativistic, so that e = cp, and that their 
fraction is small compared to the plasma density in the sheet. The singularity 
appears in the above solution of the MHD equations during the compression 
of the current sheet. The distribution function F^:(t, z,Pa:,Py,Pz) of relativistic 
particles in Cartesian coordinates satisfies coUisionless Lioville's equation 

dF, , c2 dF, , ^dF, ^ 
at s oz opi 

Let us change the Cartesian coordinates {Px,Py,Pz) to the cylindrical ones 
{p\\,p±,(p) in the momentum space. Here p\\ is a component of the particle 
momentum along the magnetic field B = B^, p± is that perpendicular to the 
magnetic field, and is a rotation angle of the particle around the magnetic 
field, 

Px = P\\ , Py =P± sin (p, pz = p± cos (f. 
The force acting on the particles is the Lorcntz force 

jr=e(E+^[pB]) . 

Thus we have the kinetic equation 

dF^ cp± dF^ ,T^. X 9F^ 

H cos (p— h e (Ey sm ^p + E^ cos ip) h 

ot p oz op± 



c 

— (Ey COS (f — Ez sin tp) 
P± 



dip 



where lOc is the cyclotron frequency of the particle rotation, lJc = eBjp. Let us 
expand the function F^, of into a Fourier series and leave only the first 3 terms 

-F* = -Fb + Fc cos ^p-\-Fs sin ip. 
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Hence we get the system of equations for the functions Fq, Fc, F^ 
dFo , cp±dF, , leEy , eEydF, , 



dt 2p dz 2 2 dpi_ 

1 eE, eE^ dF^ 

SF, cpi dFo ^ J, dFo 

OT p OZ Op± 

dFs dFo 

We need only know the function Fo{t^ z,p\\^p±), which is independent of ip, 
because the other components of the distribution function make no contribution 
after the averaging over the rotation angle (p. We solve the second and the third 
equations for Fc, Fs, and substitute the result into the first equation. This 
procedure is simplified in two cases: (a) oJcFi ^ dFj/dt, (b) UcFi ^ dFj/dt. 

In the first case, the cyclotron period is much less than the time scale of 
evolution of the distribution function. This condition is only violated close to 
the center of the sheet, where the magnetic field is zero and particles are not 
magnetized. The magnetic field is very strong in the case of the interaction be- 
tween the supernova shock and the magnetized star {B ^ 10^ Gauss), therefore 
the region, where the particles are not magnetized, is negligibly small. Hence 
we may assume that the condition (a) holds everywhere. Substituting 

F =-El dFh^cp^m p _Ey dFo 
^ B dp± eB 9z ' B dp± 

into the first equation, we obtain 

dFo EydFo cp^ d (Ey \ dFo 



— c 



dt B dz 2 dz \ B J dpA^ 



= 0. (17) 



Equation (|17p has a simple physical meaning. All particles move in the z direc- 
tion with a velocity of the electric drift —cEy/B, moreover the particles have 
an acceleration in this direction due to the variation of this velocity along the z 
axis, and hence gain the transverse momentum p±. It is important to note that 
the acceleration d{cEy/B)/dz tends to infinity at the shock point, 

d fEy \ _ Akq. ^ 



. , „ , , (l + 2fc(l-fc)g) . 

dz \ B J t 

The dependence of Ey/B on the coordinate z and the time t is given by the 
formula ()16p . Now, it will be convenient to change variables from {t,z,p±) to 
{rj,p,p±). In the new variables p7|) reads 

(fc + 2),— -p— --p^— = 0. (18) 

We may assume that the distribution function of fast plasma particles is Maxwellian 
at t = (77 = 1), and their density is proportional to the plasma density, i.e. the 
fraction of the accelerated injected particles is constant along the entire sheet 

Fo{v = l,p,P±)^cope~'P^^^". 
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The solution of kinetic equation (jlSp with such initial conditions is 
i^o(^,P,P±) = copr;i/('=+2) exp (-cp^^V(2fe+4)/y^^ _ 

It shows that the distribution function of the particles remains Maxwellian, but 
the temperature is time-dependent T — Torj^^^'^^'^'^'^^ For —2 < fc < the 
temperature is rising as rj tends to zero. At the moment 77 = the tempera- 
ture becomes infinite, which means that the distribution function is constant in 
the region p± < 00. Note that the moment j] — occurs before the shock for 
— 1/2 < A: < and these moments coincide if fc = ^1/2- The case k — —1/2 cor- 
responds to the initial distribution of density p{z) with constant initial plasma 
temperature along the entire sheet. Though formally the particles are acceler- 
ated up to infinite energies, it is clear that the cyclotron radius of the fastest par- 
ticles should be less than the typical thickness of the sheet L, and also, that their 
cyclotron frequency should be greater than the inverse typical evolution time At 
of the sheet. Therefore, we have two estimates for p±rn- P±m = Pim = eL < B > 
I c and cp±m — cp2m — ToATy""'^/^^'^"'"'*^ , where Arj can be found from the defini- 
tion of 77 and is estimated to be A77 ~ {cs/{Lujco))'^{p2m/'n^c)'^, = eB/mc. 
Finally we have P2m/mc = (To/(toc2))('=+2)/('^+3) (Lt^^o/cs)^/'''+^^ ■ Later we 
will see that pim >> P2m for the actual values of the parameters, so we can 
assume that p,„ = p2m ■ 

5 Synchrotron radiation 

The most interesting case for us is when fc ~ —1/2, and in this case all plasma 
masses up at the moment of shock near the two planes = ±2/oy^, 

symmetric with respect to the center of the sheet. The estimate for z+ is 
2Cso/l^''^/^''llt=Oi it is reciprocal to the initial acceleration of the sheet dv/dt\t=Q, 
and, of course, z+ should be less than the thickness of the sheet L. The plasma 
velocity Vz is positive for z > and Vz < for z < 0, i.e. plasma particles move 
away from the center of the sheet into the regions around z+,- with a strong 
magnetic field. Due to the strong magnetic field, the synchrotron radiation be- 
comes significant for the plasma particles (electrons) there. For z > z+ and 
z < z_ the quasiclassical approximation used in the solution of the MHD equa- 
tions is not applicable. We assume that the magnetic field is almost constant 
at the large distances z > z+ and z < z_ from the center of the sheet, and its 
value is equal to the magnetic field at the boundary of the sheet. The electric 
field with the components Ey and Ez exists there as well, and plasma velocity in 
the electromagnetic field should be w ~ Cs, where Cs is the fast magnetic sound 
speed. 

Thus, fast particles with the constant distribution function obtained above 
move into the regions with constant fields due to the electric drift and strongly 
radiate there. Next, we find F, the distribution function of these particles, and 
their radiation intensity /. 

We may assume that Ey/B ^ 1 and Ez/B ^ 1, i.e. the plasma motion is 
non-relativistic. In this case the drift velocity is v^^ — (0, cEz/B, —cEy/B). 

The function F satisfies the following continuity equation in the phase space 

dF _^ ^pz dF ^dTi^^^dF _^ 
dt p dz dpi ^ dpi 
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The force Ti acting on the particles has the simplest expression in the reference 
frame moving with the velocity of the electric drift, where there is no electric 
field. Since the drift is non-relativistic, this force is the same in the stationary 
reference frame. We use tildes for the values observed in the moving frame. The 
friction force due to radiation is [6] 

= -^^l (pI+pI) iP-^Pv^Pz) , 

note that the magnetic field in the moving reference frame is i? = i? (l — [Ey + El) / B^) , 
and components of the electric field are Ey = 0, E^ = 0. Let us denote 

2e4 



The net force acting on the particles is 

^ X -F fr^x •} 

*Fy cB—^ -\- J~ fr^yi 

J'z = -eB^+Tfr,z- 
P 

Again, change the Cartesian coordinates {px,Py,Pz) to the cylindrical ones 
in the momentum space 



p± sm if 



Px 


= Px ^ P\\ 






Ez . . 


Ez 


Py 


= ^P + Py = 


-bP + p 








Pz 


--^P + P^ 


= -p 

B ^ 



Notice that a non-relativistic momentum addition rule was used here as the drift 
is non-relativistic. 

Further, we neglect terms of the order higher than Ey/B and Ez/B in all 
calculations. Under this approximation B = B, 

^ OF _ ^ dF TyPyEz + TzPzEy dF ^ J'yPyEy + TzPzEy dF ^^g^ 
^ dpi 'dpi Bp dpy Bp dpz 

Since divergence is a scalar value and the Lorentz force is solenoidal, we have 

dTi _ dTi _ dTfr,t _ _4^^2Pi 



dpi dpi dpi p 

Substituting the expressions obtained above in the original kinetic equation, 
and averaging it over the rotation angle Lp, we obtain 

dF cEydF ,p2 / OF \ p2 

dt B dz P \ dp\\ dp± ) p 

The identity p = p is used here, it can be proved by expanding p into a Taylor 
series in powers of Ey/B and E^/ B and averaging over the angle (p. 
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Let us consider the case z > (z < is analogous due to the symmetry 
of the problem). The particles with constant initial distribution function move 
into the region z > z+ after accumulating at z+ = 2|/o|y^, according to the 
kinetic equation obtained above. We can suppose without loss of generality that 
the origin of the z coordinate is at the point z+, i.e. we can solve the kinetic 
equation with the initial condition 

F(t = i»; z = 0;p||;p±) = cie(p„ -p). 

Here is the time of the constant distribution formation. Let us notice that 
that the constant distribution function of accelerated particles obtained in the 
previous section was the function of the transverse momentum Here, in 
contrast, we use the isotropic function with the same parameters. The reason 
is that the plasma with strongly anisotropic particle distribution is very un- 
stable with respect to the excitation of the electrostatic waves leading to fast 
isotropization. The growth rate of such instabilities is of order of the cyclotron 
frequency u)c or the plasma frequency Wp jlOj . We are looking for a stationary 
distribution function (the solution of ([20|) which does not depend on time) to 
find the spectrum of the radiation of the particles. This solution is 

F{z]p\\]pi_) = , e pm 



cEyP 



Having found the distribution function, we can compute the spectrum and the 
directivity of the radiation, i.e. dependence of the radiation intensity on the 
direction and on the frequency. 

First, let us find the synchrotron radiation from one particle moving with the 
velocity v. The prime denotes that the variable is observed in the coordinate 
frame K' moving with the velocity U|| . To find the radiation, we will need 
to change the coordinate frame from the stationary one to K' and back. We 
denote the angle between the direction of the magnetic field and the velocity of 
the particle by x and the angle between the plane of the particle rotation and 
the direction of the radiation by 9. The radiation in the K' coordinate frame is 



2 - cos^ 9' - ^ [1 + ^] cos* 9' 



Here 



1 1 — 73- COS"' 9' j 
dVt' = 2TTCOs9'd9', B' ^ 7S, v'j^ - 7wj_; 



v±_ = u sin X, w|| = w cos X, 7 = 1/ V^~(3^, P c. 
Therefore, the radiation in the stationary coordinate frame is 

dl ^ {1 + Psin9')dr, 

. [3 + sine' 

sm6^= , (21) 

l + /5sin6l'' ^ ' 

d9 , ^ d9' 

- 7 l + /5sin0O^757' 
cos-^ 9 cos-^ 9' 
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(COS 



dJ _ ^ 

d6 ^ cos^ ( 



3«' e^B^v l{l~-i''vi/c') 



2 - cos2 e 



2 2 / 

4c2 



COS^ 9' 



7/2 



1 -72^C0S2 fl'' 

The sought dependence of the radiation intensity on the direction is 



^de' 



dz 



p dp 



1 



cos^ 6 



2 - cos^ 


2 2/ o 2 2 \ 


cos^ 61' 




(l-72^COs2 0') 


7/2 



Cl 



^psin^x 



1 + ^Psin^x 



X 2iT smxdx- 



(22) 



Note, that only the particles with 7 3> 1 make significant contribution to the 
radiation. In the K' coordinate frame the synchrotron radiation is concentrated 
near the plane of the particle rotation, which means that 9' = Q. The most part 
of the radiation is concentrated within the angle /S.9' ~ 1/7^ 1. Note that 
r = (1 - i;2/c2)-i/2 ^ 1 since 7 > 1. 

From Equation (|2ip it follows that sin6' ^ (3 = ^cosx ~ cosx. Hence 
the particles with the pitch angle x radiate in the stationary coordinate frame 
mainly into the angle 

(23) 

Let us perform the Taylor expansion of all expressions appearing in the integral 
over X in (|22p in powers of 1 — d^/c^ and leave only the first, most essential 
terms. Using the smallness of the angle 9' , we get 

7 ~ 1/sinx, cos 6*' ~ 1, ~ csinx, 
2 X 1 



i ^ sm X 
7 



? 2 

.2 r^i 

^ " 7^ 

4c2 



^]cos^.' 



4p2 „;-,2 



sm X 



Substituting these equations into the integral (|22p and using ([23]) . we obtain 
that the integrand is approximately equal to 



7 e^B^ 
77C1 



8 m^c'^ 



1 + ^^^p cos2 ( 



X e p„ 



kB±z 



p cos2 6* 



dx- 



Numerical computation shows that the radiation at the angle 9 is produced 
mainly by the particles with Ax — 1.52/r, i.e. the integral over dx is close to 



1.33ci- 



cos 9 



^pcos^i 



X e p,n - 



1 



p cos^ 9 



13 



The upper limit of the integration over p in Equation 



is 



P -Per - ^g3^ 



l-^PraCOS-'e 

Notice that pcr < Pm since Ey/B < 0. The lower limit of the integration, 
can be replaced by zero. Hence 



— ) = 1.33ci — ;r^cos6' / dz P x IH — pcos 6* dp. 

d0 m'^& Jq Jo V '^Ey 

These two integrals can be evaluated analytically: 

(-) = 1.33c, f^)%,_^ (3K^^ 
^dO' Um'^c^ \ K ) cose ^ 

4K^ 6K^ UK 12 



cos4 e cos6 9 cos8 e 



log(l + Xcos^ 0) 



where K = —(2e B /Sm c Ey)pmZmax 2> 1. Angular distribution of the inten- 
sity {dl /dff) is maximal at cos 6* ~ K^^^"^, it vanishes at cosfl = 0, and has the 
width maximum A(cos0) close to K~^/^. That is, the radiation is strongly col- 
limated along the magnetic field and the coUimation angle can be estimated 
as K~^/'^ <^ 1. Let us find the radiation spectrum. For a single particle the 
spectral intensity of the synchrotron radiation is [7] 



dl V3 e^Bj 



tjj 



F — 



duj 2ti mc^ V ujc 

Here 



P smx, 



2mc \mc^ J 2m?c^ 
B± — Bsinx- 

The total spectrum is given by 

,dr . , /■^"- , /-p- / kB'^z .2 V^di o, 

(— )=ci / 2'Ksmxdx dz X / 1 + — — psm x -i-P dp, 

duj Jq Jq Jo \ cEy J duj 



(24) 



where 



kB'^z 2 
— — p™sm X 

CEy 



Asymptotic behavior of the function F{^) is known to be 



m-\l^^-' ase 

Using Laplace's method, we have the asymptotics when uj ^ oo 
/ K.B^z 2 dl 3 

e'^B^ V^^ 4 ./2 / 2m^c^ 1 

p„sm ' xexp -uj- 



3eB p^^sinx 
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After substituting the expression for pcr, we obtain the asymptotics for the two 
other integrals when cj — > cx) 



eB'' sm X j \ 6eB sm x 

Finally, we obtain the asymptotic of the spectral intensity 



,dl 277r\/3 6 e\Ey\ 1 / 



The exponential damping of the spectrum begins at the frequency 



3 eB /PmY 
2 mc \mcJ 



Now let us find the asymptotics of the spectrum as a; ^ 0. Recall that the radi- 
ation is mainly concentrated in the angle K~^^^ radian. Only the particles with 
X ^ K^^l"^ make a substantial contribution (equation ([23]) ). The asymptotics 
of F{i) is 

Therefore, when lj ^ 0, we have 



Finally, after integrating over x. and z, we obtain 

27 ^ ^ /7\ m4c6|i;j,| 4,3 1/3 



'duo' ~ '14 V6/ e4/3S7/3 

The asymptotic behavior of (dl/dcj) changes from lo^^^ to uj^^ somewhere in 
the region w < . An approximate value of the frequency uj ~ loq when this 
change occurs is given by 

mc Kmc/ 

Thus we see that the radiation spectrum of the accelerated particles of the whole 
sheet is maximal at a; ~ cjg with power gain for lo < uj^ and power damping for 
Lo > uj^. The exponential damping occurs for the large frequencies > w*. 



6 Summary 

According to the papers of the typical moderate values for the sheet are: 
the typical velocity of the fast magnetic sound is ~ 10^ cm/s, 
the thickness of the sheet is L ~ 10^ cm, 
the typical value of the magnetic field is B ~ 10^ G, 
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Figure 3: Radiation spectrum 



the plasma temperature in the sheet is Tq ~ 10 keV. 

The time of the sheet compression can be estimated from Equation (jisp and 
the condition z+ ^ L 



a 



dv 
dt 



L 

— ^ O.lsec. 



Now we can estimate the maximum value of the Lorentz factor for the particles 
(fc = -l/2): 

/ \ 3/5 / r \ 2/5 



mc 



10" 



The collimation A9 of the radiation is determined by the parameter K, A9 

2 Trnio'^oLre B 



K 



3 C2 \Ey\ 



Tf. is the classical electron radius, rg = 3 • 10^^'' cm. Assuming \Ey\lB ~ Cg/c, 
we obtain K ~ lO'', which yields the collimation angle, AQ ~ 3 • 10^^. The 
intensity of the accelerated particle radiation in the sheet is maximal at the 
frequency cjq — 8 • 10^*^5^^, which corresponds to the energies of the order of 
500 keV. Asymptotic behavior of the spectrum as well as the frequency agree 
with the observed values of GRBs [2j. 

Thus, we see that the more precise calculations for the model of cosmolog- 
ical GRBs proposed by [S] confirm their estimations. Evolution of the current 
sheet, particle acceleration due to the magnetic energy release, and the syn- 
chrotron radiation turn out to agree with what was estimated. The total ra- 
diated energy corresponds to the release of the magnetic energy stored in the 
sheet (10'^'' — lO^'^erg, see the papers of [S]). The equivalent radiated energy in 
the isotropic models is (3 • 10^)^ times greater, i.e. about lO''^ — 10'*''erg. The 
latter energy can be even greater for closer binary systems (where the distance 
between a supernova precursor and the star-companion is less than 10^"^ cm), 
and it may reach values of lO'^" — lO^^erg. The proposed model explains the 
main characteristics of observed GRBs. 
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